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The Jackiw-Rebbi model describes a one-dimensional Dirac field coupled to a soliton field and can be 
equivalently thought of as a model describing a Dirac field with a spatially dependent mass term. Neglecting 
the dynamics of the soliton field, a kink in the background soliton profile yields a topologically protected 
zero-energy mode for the field, which in turn leads to charge fractionalisation. We show here that the model, 
in the first quantised form, can be realised in a driven slow-light setup, where photons mimic the Dirac field 
and the soliton field can be implemented-and tuned-by adjusting optical parameters such as the 
atom-photon detuning. Furthermore, we discuss how the existence of the zero-mode and its topological 
stability can be probed naturally by studying the transmission spectrum. We conclude by analysing the 
robustness of our approach against possible experimental errors in engineering the Jackiw-Rebbi 
Hamiltonian in this optical setup. 

Quantum simulators first and foremost offer a promising alternative when analytical and numerical 
methods fail in analysing models with strong correlations. On the other hand, they can also be used 
in probing exotic physical phenomena such as those predicted by relativistic theories. To date, a col- 
lection of effects in different fields ranging from condensed matter physics to relativistic quantum theories and 
material science have been simulated\ using different platforms such ion traps^, and cold atoms in optical 
lattices^. 

More recently, the ability to controllably manipulate photons and their interactions with atomic systems 
resulted in the birth of a new direction in quantum simulations using photons and polaritons to mimic strongly 
correlated phenomena. Coupled cavity QED arrays (CCAs) were initially considered, where photons trapped in 
resonators interfaced with two level atoms (real or artificial ones) were shown to be able to reproduce many-body 
dynamics^"^. The so-called photon blockade effect was exploited and polariton Mott transitions and effective 
spin-models were proposed, introducing what is now known as the Jaynes-Cummings- Hubbard (JCH) modeF. 
Simultaneously and independently, the possibility for strongly correlated polariton dynamics in CCAs with four- 
level atoms and external fields was proposed^, followed soon after by the JCH's phase diagram^'^° and the 
Fractional Hall effect^ \ More recently one dimensional highly nonlinear waveguides with slow- light nonlinea- 
rities^^"^'* have been considered where effects characterising Tonks- Girardeau gas, Luttinger liquids physics, or 
even interacting relativistic theories have been shown to be simulable^^"^^. The possibility to probe out-of- 
equilibrium phenomena has also been explored in driven setups^^'^^ 

Quantum simulation has also shown great development in bringing the exotic physics of single particle 
relativistic effects into laboratory. These range from numerous theoretical works studying the Dirac equation 
and emerging effects such as the trembling motion of the electron (Zitterbewegung) or Klein tunneling to 
experimental implementations of those in different platforms covering all three mentioned above. More specif- 
ically in ion technologies the Dirac equation in 1 + 1 dimensions has been realised^""^^ followed by an imple- 
mentation with photons in waveguide arrays, including the random mass Dirac modeP^"^^. Seminal proposals for 
the realisation in slow- light systems also exist^^"^°. Going beyond the Dirac equation, recent works have proposed 
to simulate the Majorana equation^^'^^ and neutrino oscillations^^ in trapped ions. 

The purpose of this article is to propose a quantum simulation of a historically important relativistic model 
known as the Jackiw-Rebbi (J-R) modeP^ with slow-light. The model describes a one-dimensional Dirac field 
coupled to a static background soliton field and can be equivalently thought of as the model describing a massless 
Dirac particle under a Lorentz scalar potential. The same model has been studied independently by Su, Schrieffer, 
and Heeger (SSH), while studying electron-phonon coupling arising in polyacetylene, in the continuum limit^^. 
The model is well known for predicting charge fractionalisation^^, well before fractional quantum Hall effect was 
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discovered, and also for the topological nature of its zero- energy 
solution and can be thought of as a precursor to topological insula- 
tors, a topic that is being hotly pursued^^'^^. 

There have been proposals to realise the model and observe charge 
fractionalisation in an optical lattice setup^^'^° and experimental 
observation of the soliton which follows the model in a fermionic 
superfluid^\ Also, there have been proposals to observe topological 
bound states analogous to that arising in the model, in quantum walk 
and graphene setups^^"^^. A topological bound state has also been 
observed in optically- induced dimer lattices^^, where the latter is 
closely related to the SSH model. Here, we follow a different route 
and propose a photonic implementation of the model in a slow-light 
polaritonic setup and show that the topological properties can be 
probed straightforwardly in an optical transmission experiment. 
We would like to highlight here that the usual experimental difficult- 
ies, in realising the photonic nonlinear interaction in slow-light sys- 
tems for more complex many body simulations, do not exist in this 
case. Therefore we believe the proposal is a good candidate for 
directly realising the J-R model and could allow the efficient probing 
of its topological properties for the first time. 

The Jackiw-Rebbi model and the related model studied by Su, 
Schrieffer, and Heeger^^ share many features similar to those studied 
in topological insulators and might indeed be classified as an AIII- 
type chiral topological Dirac insulator under a suitable regularisa- 
tion^^. Here we describe the Jackiw-Rebbi model and point out the 
topological properties and similarities to topological insulators. The 
equation of motion considered by Jackiw and Rebbi reads 



idt^ = dcpz + 



^(z) 



1 



where a and are the Dirac matrices which in this case can be chosen 
to be proportional to two PauU matrices. For concreteness, we take 
them to be a = Gz and = Oy. The real scalar field <l){z) is assumed to 
obey the Klein-Gordon equation of the form 



2k;2 



(2) 



The ground state of the scalar field is degenerate at (l){z) = ± k, which 
implies the existence of a soliton that interpolates between —KdXz = 
—00 and k: at z = 00 and corresponding anti-soliton. The soliton 
localised at z = 0 is described by 



^{z) = K tanh(>^2:) 



(3) 



When the Dirac field is coupled to such a soliton, a zero-mode (zero 
energy state) appears which is localised around the soliton. The 
unnormalised spinor wave function of this zero-mode reads 



1*0(2) = exp 

= exp 



dx(l)^ [x) 



1 



(4) 



■ln(cosh kz) X- 



where a^x ^ ~^X^ ^^d is shown in Fig. 1 (soUd red curve). Note that 
the dynamics of the scalar field is neglected in the above argument, 
i.e., the scalar field is treated as a constant background field. 

Far away from the kink, the particle and hole bands have an energy 
gap m, whereas this gap must close at the point (l)s{z) = 0. This 
resembles the gap closing at the boundary of a topological insulator 
where the bound surface mode develops. The resulting bound state, 
the zero -mode, is protected by the topology of the scalar field, whose 
existence, irrespective of the local profile of the kink, is guaranteed by 
the so-called index theorem^^'^^. This phenomenon is similar to the 
emergence of edge modes in the quantum Hall effect^^'^° or topo- 
logical insulators^^'^^, where edge modes appear at the boundary of 
two topologically different domains. 
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Figure 1 | Soliton profile, ^^(z) (blue dotted), and the zero-mode wave 
function profile, |^o(^) |^ (red solid), showing the localisation scale A of 
the zero-mode. We have set m = c = /c = 1 for convenience. 

Another interesting aspect of the model (when second quantisa- 
tion is taken into account) is charge fractionalisation, which we 
briefly describe before moving on to a proposal for a photonic imple- 
mentation. The ground state of the Dirac field (the vacuum) in the 
soliton background may or may not include the zero -mode. Because 
of the charge conjugation symmetry of the system, the two degen- 
erate ground states must have opposite charges. Moreover, the 
charge difference between the two must be 1, as there can be only 
one electron occupying the state. The result is that the filled or 
unfilled states must have charges 1/2 and —1/2, respectively. This 
result can be confirmed by constructing the formal charge operator 
in terms of creation and annihilation operators of the eigenmodes. 
The vacuum states are the eigenstates of the charge operators which 
means that the observed fractional charges are sharp and are not just 
a trivial realisation of a distributed charge. 

Results 

Slow-light realisation of J-R dynamics. The spinor slow- light setup 
we employ to realise the J-R model and subsequently observe the 
zero-mode is depicted in Fig. 2. This setup was first proposed by 
Ruseckas et al.^°, to which we refer the reader for a detailed 
explanation. The system comprises of a waveguide system coupled 
to an ensemble of atoms, which could be realised either in a tapered 
fiber^^'^^ or a hollow-core waveguide^^"^^. The atoms are characterised 
by the three hyperfine ground levels; one populated state |^) and two 
unpopulated states |sl) and |s2). These ground states are coupled to 
the two excited states \el) and \e2) by the probe and control fields. 



Si{z = 0) = a 



Si{z = L) =aT 




E2{z = 0) = aR 




Figure 2 | (a), Schematic diagram of the optical waveguide system, 
interfaced with an ensemble of atoms where propagating light fields Si and 
S2 play the role of Dirac spinor components. By adjusting the relevant 
optical couplings and detunings, the J-R model can be simulated and its 
topological aspects probed by looking at the transmission spectrum, 
(b), The level structure of the interfaced atoms. 
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The counter-propagating probe beams are described by the electric 
field amplitudes Ei and E2, with the respective central frequencies coi 
and CO2 and drive the transitions |^) ^ \el) and |^) \e2). The 
propagation of the probe beams is controlled by two pairs of 
counter-propagating control lasers with Rabi frequencies Qji and 
Qj2 (with j = 1, 2), driving the transitions from the excited states 
to the unpopulated states. A slowly varying amplitude Sj{r,t) is 
associated with the electric field strength Ej of thej^^ probe field: 



Ej{z,t)- 



lhco~, . 

-—Sj{z,t)e 
^0 



ht+ikj2 



-\-c.c., 



(5) 



with ki = Wi and k2 = —(JO2, where the speed of light in an empty 
waveguide is taken to be 1. 

The propagation of the slowly varying amplitudes is such that they 
follow the 1 + 1 dimensional Dirac-like equation 



1 + 



1 1 
Vo sin^S 
S 



1 cos S 
Vo sin^S 



dtS + d^S 



(6) 



Vo sin S 



Vq = Q^/g^n is the group velocity much smaller than 1, where Q j \fl 

is the Rabi frequency of the control fields,^ is the atom-light coupling 
strength, and n is the atomic density. The complex Rabi frequencies 
are tuned so that Sn = S22 = 0 and S12 = ^21 = 5 where 

Qij = Q j \/2exp(/S2)). ^is a column vector of two slowly- varying 

field components E\ and E2, i.e., E= (^1,^2)^ and b is the two- 
photon detuning. In the limit S = 7i/2, the above equation reduces 
to the Dirac equation in 1 + 1 dimension 



[idt + ivoGzdz 



i7) 



which is equal to the Dirac equation introduced in the previous 
section with the identifications c = Vq and mc^lK = 3. 

The connection with the J-R model is obvious once we let 3 be 
spatially varying as 3{z) = 3o tanh(/lz). This is possible, for example, 
by introducing additional lasers and states to induce AC Stark shift in 
the metastable states |si) and |s2). Spatial variations in the lasers can 
then induce an arbitrary mass profile, including the desired soliton 
profile. A similar model has been studied in Ref 29 with a slightly 
different atomic level scheme and the existence of zero-mode has 
been briefly commented on, although no connection with the J-R 
model has been made nor the topological nature of the zero -mode 
mentioned. By making the connection, it is easily seen that there is 
interesting physics to be explored in the slow-light system, namely 
the topologically protected zero-mode. In the following, we discuss 
how the zero -mode and its topological stability can be observed 
experimentally. 

Probing the zero-mode and its topological stability. Broadly 
speaking, there are two possible ways to observe the zero-mode in 
this optical setup. The first is an adiabatic method, where a photonic 
wave packet is adiabatically loaded to prepare an initial state that 
resembles the zero-mode whose evolution is then observed. In this 
method, the initial state is prepared by capturing a pulse in the 
medium via the usual electromagnetically induced transparency 
(EIT) technique, i.e., by slowly turning off the forward-traveling 
control fields which are initially on^^. Then, all the control fields 
are slowly turned on again, including the coupling fields Q12 and 
^21- At this point, the pulses are trapped as stationary light and go 
through the effective evolution governed by the Dirac equation. The 
dynamics of the spinor fields then differ significantly, depending on 
the initial condition of the wave packets as shown in Fig. 3, where we 
have solved the Dirac equation with the initial gaussian wave packet 




Figure 3 | Dynamics of an optical gaussian wave packet mimicking the 
zero-mode wave function (4) in the J-R model (top row), compared with 
its evolution under the Dirac equation (bottom row). The overlap between 
the zero-mode and the gaussian wave functions is ~ 0.997. In the first case 
the coupling to the background soliton forces the initial wave packet to be 
trapped, while the second case shows the expected wave function spreading 
for a Dirac particle of mass ^o- ^1 ^i^d S2 corresponds to the top and bottom 
parts of the spinor, respectively. 



exp[ — z^/2o-^] ^ Y V^o-]/ with 0- = 1.2 and/ = (1, 1)^ to mimic the 

zero-mode spinor. In the case of the J-R model, i.e., under the soliton 
mass profile, the wave packet is trapped around the kink (upper row 
of Fig. 3), while under the normal case of constant mass, the wave 
packet spreads out (bottom row of Fig. 3). 

The second method, better suited for this particular experimental 
realisation, is to look at the transmission and reflection of an incident 
probe field Si{z = 0). Similar studies have been carried out quite 
recently, using transmission to probe strongly interacting effects in 
similar polaritonic systems^^"^". Consider a monochromatic probe 
field Si = ae~^^^^ impinging from the left, while S2{z = L) =0 where 
L is the length of the system. One can study the transmission and 
reflection spectrum of this field where the transmitted field will come 
out to the right of the waveguide as Si{z = L) = aTe~^^^\ whereas 
the reflected part is S2{z = 0) = (xRe~^^^^ as shown in Fig. 2a. The 
constant mass case can be solved analytically^^ and shows the beha- 
viour depicted in Fig. 4a. There is a finite window of perfect reflection 
due to the well known band gap proportional to the mass energy. 
Upon introducing the spatially dependent mass term discussed earl- 
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Figure 4 | The reflection \R\^ (black) and transmission | (red) curves 
for the effective Dirac particle (a), without the soliton background and 
(b), with a soliton field whose profile is 0.25 tanh(0.02z). (a) shows the 
Dirac mass bandgap whereas (b) shows near-unity transmission near the 
zero-energy due to the bound zero-mode. 
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Figure 5 | (a), (b), Stability of the zero-mode under changes in the background profile, and (c), (d), random fluctuations in the soliton profile, 
implemented via the two photon detunings. (a). Background profile changed to the sine function d{z) = 0.25 sin(O.Olz). (b), Soliton profile changed to 
the linear fijnction going from —0.25 to 0.25. In both cases, the response of the system around resonance remains similar to Fig. 4b. (c), (d). Added 
random fluctuations where the maximum value of the random fluctuation e(z) is 0.2 and 0.5. We see that the zero-mode persists in spite of the relatively 
large values of fluctuation although the maximum possible transmission is gradually reduced from roughly 95% to 85%. 



ier, i.e. tanh(>^2:), transmission reappears at the center of the bandgap 
as shown in Fig. 4b, which can be attributed to the presence of the 
zero-mode. Here and below, we foflow Ref. 30 and assume that the 
length of the atomic cloud interacting with the propagating light is L 
= 300 iim and the latter's group velocity Vq = 17 mis. For these 
values, the maximum two photon detuning — 0.25vo/L lies wefl 
within the EIT transparency window. 

Using this second method, the topological nature of the zero- 
mode can be readily probed by looking at the transmission spectrum 
whfle perturbing the mass profile. The latter can be done by tuning 
the two photon detunings using standard optical methods, like AC 
stark shifting used in slow-light experiments Changing the 
hyperbolic tangent function to the sine function whfle preserving 
the topology of the profile, i.e., the mass term takes the value of 
— ^0 and +dQdX z = —LI 2 and z = LI 2 respectively and crosses 0 
only once at the origin, has little effect on the transmission at the 
center of the bandgap as shown in Fig. 5a. The extreme limit where 
the 'kink' has been straightened out such that the background profile 
is linear is also shown in Fig. 5b. The zero-mode persists in this 
extreme case, clearly demonstrating the topologically protected 
robustness. The zero-mode is also protected from random fiuctua- 
tions in the mass proffle as shown in Fig. 5c and d, where we have 
assumed random fiuctuations within 20% and 50% of d{z), i.e., 
6{z) = ^otanh(0.022:)(l +e(2:)) where e(z) is a uniform random with 
the resolution of 0. 1 ^m. Such topological protection holds as long as 
the fiuctuations do not destroy the topology (two distinct values at 
the infinity and one zero- crossing only) of the soliton background. 
Unusually large fiuctuations that changes the topology could destroy 
the zero-mode, but are unlikely in experimental implementations. 

Experimentally, it might be difficult to set the phases of the control 
fields exactly to the required value of the mixing angle S = nl2 in Eq. 
(6), in which case the resulting equation deviates from the Dirac 
equation. To study the effect of this imperfection, we have set S = 
7i/2 + AS where AS takes 10% and 40% of the desired value, the result 
of which is depicted in Fig. 6a and b, respectively. The presence of the 



zero -mode persists upon experimental errors in creating the exact 
Hamfltonian. Another experimental imperfection that can affect the 
transmission spectrum is the loss. The loss from the excited states can 
be easily incorporated into the equation of motion by introducing an 
extra term — /yg^into Eq. (7)^°, where y^ff = yS^IQ^ with y the decay 
rate from the excited states. Effects of photon losses are shown in 
Fig. 6c and d, where we have used 3(z) = Sq tanh [0.022:] and 
y^Q/Q^ = 0.005 and 0.01, respectively. The transmission peak at 
the zero -mode reduces in size with increasing loss rate, but the peak 
is stifl visible at ydl /Q^ = 0.01. Under normal EIT condition, 3o<^Q, 
which along with our numerical analysis indicates that the effect of 
the loss should not obstruct the detection of the zero-mode. 
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Figure 6 | (a), (b). Departure from the Dirac equation for different values 
ofthe mixing angle (seeEq. (6)): AS = 0.l7i/2 and0.47i/2. (c), (d). Effects of 
loss from the excited states for y^g/Q^ =0.005 and 0.01. The zero-mode 
peak survives all these experimental imperfections. 
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Discussion 

In this work, we have shown that it is possible to simulate the Jackiw- 
Rebbi model and probe its topological nature in a driven slow-light 
system using current technology. By introducing spatially- 
dependent optical detunings, it is possible to simulate the Dirac 
equation with a spatially dependent static soliton field, and allows 
one to directly probe the topologically protected zero-mode. The 
robustness of this zero -mode can be tested by changing the spatial 
profile of the detunings, while continuously observing the optical 
transmission spectrum of the system. We have also analysed the 
effects of experimental imperfections such as deviations from the 
ideal model and losses and showed that the observation of zero -mode 
survives these imperfections. 

Before closing we would like to briefly comment on other types of 
models that might be studied in the same system in future works. 
Firstly, adding interactions in the system through EIT photon non- 
linearities could allow the study of interesting strongly correlated 
physics. Intra-species repulsion would make the bosons behave like 
fermions in many ways, and it would be interesting to think about 
how this affects the (lack of) charge fractionalisation effect in bosons. 
Interacting random mass Dirac model is another interesting possibil- 
ity as studying the interplay between interactions and randomness is 
an important field being actively studied. Yet another possibility is to 
think of the mass term as the Lor entz- scalar potential which can act 
as a confining potential, given a proper spatial dependence^^'^^. This 
model has been shown to act as a phenomenological model of quark 
confinement motivated by the 'MIT-bag' modeP^. 
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